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Abstract: It is shown that the dissipation due to spontaneous emission can entangle two closely 
separated two-level atoms. 
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1 Introduction 



Analysis of various aspects of spontaneous emission by a system of two-level atoms, initiated in 
the classical paper of Dicke pP was further developed by several authors (see e.g. [21 El El)- In 
particular, in the case of spontaneous emission by two trapped atoms separated by a distance small 
compared to the radiation wavelength, where is a substantial probability that a photon emitted by 
one atom will be absorbed by the other, there are states of the system in which photon exchange 
can enhance or diminish spontaneous decay rates. The states with enhanced decay rate are called 
superradiant and analogously states with diminished decay rate are called subradiant pQ. It was 
also shown by Dicke, that the system of two coupled two-level atoms can be treated as a single four- 
level system with modified decay rates. Note also that such model can be realized in a laboratory by 
two laser-cooled trapped ions, where the observation of superradiance and subradiance is possible 

m 

Another aspects of the model of the spontaneous emission are studied in the present paper. 
When the compound system of two atoms is in an entangled state, the irreversible process of 
radiative decay usually destroys correlations and the state becomes unentangled. In the model 
studied here, the photon exchange produces correlations between atoms which can partially over- 
come decoherence caused by spontaneous radiation. As a result, some amount of entanglement 
can survive, and moreover there is a possibility that this process can entangle separable states of 
two atoms. The idea that dissipation can create entanglement in physical systems, was recently 
developed in several papers El EJ ■ In the present paper we show that the dissipation due to 
spontaneous emission can entangle two atoms that are initially prepared in a separable state. We 
study the dynamics of this process. In the Markovian approximation it is given by the semi-group 
{T t } of completely positive linear mappings acting on density matrices ^U]. We consider time 
evolution of initial state of the system as well as the evolution of its entanglement, measured by so 
called concurrence El j m the case when the photon exchange rate 7 is close to spontaneous 
emission rate of the single atom 70 and we can use the approximation 70 = 7 (similar model was 
also considered in ^3]). We calculate asymptotic stationary states p as for the semi-group {T t } and 
show that they depend on initial conditions (i.e. {T t } is relaxing but not uniquely relaxing). The 
concurrence of p as also depends on initial state and can be non zero for some of them. We discuss 
in details some classes of initial states. In particular, we show that there are pure separable states 
evolving to entangled mixed states and such which remain separable during evolution. The first 
class contains physically interesting initial state when one atom in in excited state and the other 
is in ground state. The relaxation process given by the semi-group {T t } produces in this case the 
states with entanglement monotonically increasing in time to the maximal value. The class of pure 
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maximally entangled initial states is also discussed. Similar "production" of entanglement is shown 
to be present for some classes of mixed states. On the other hand, when the photon exchange rate 
7 is smaller then 70, the relaxation process brings all initial states to the unique asymptotic state 
when both atoms are in its ground states. As we show, even in that case the dynamics can entangle 
two separable states, but the amount of entanglement is decreasing to zero. 



2 Pair of two- level atoms 

Consider two-level atom A with ground state |0) and excited state This quantum system can 
be described in terms of the Hilbert space Ha = C 2 and the algebra 21a of 2 x 2 complex matrices. 
If we identify |1) and |0) with vectors (J) and (J) respectively, then the raising and lowering 
operators er + , cr_ defined by 

a + = |l)(0|, a_ = |0)(l| (1) 
can be expressed in terms of Pauli matrices o~\ , o~i 

a + = ^(a 1 +ia 2 ), a- = ^ (01 - i ct 2 ) (2) 

For a joint system AB of two two-level atoms A and B, the algebra %ab is equal to 4 x 4 complex 
matrices and the Hilbert space Hab — Ha <& Hb = C 4 . Let £ab be the set of all states of the 
compound system i.e. 

S ab = {p e %ab : p > and trp = l} (3) 
The state p S Sab is separable |14| . if it has the form 

p = J2 Xk Pk®Pk> p£^£a, p% e£ B , Xk>0 and J2 Xk = 1 W 

k k 

The set £^b 01 a ^ separable states forms a convex subset of Sab- When p is not separable, it is 
called inseparable or entangled. Thus 

£ab = £ab \ Sab ( 5 ) 

If P £ Sab is a pure state i.e. P is one-dimensional projector, then P is separable iff partial 
traces tr^P and tr^P are also projectors. For mixed states, the separability problem is much 
more involved. Fortunately, in the case of 4 - level compound system there is a simple necessary 
and sufficient condition for separability: p is separable iff its partial transposition p TA is also a 
state |15j . Another interesting question is how to measure the amount of entanglement a given 
quantum state contains. For a pure state P, the entropy of entanglement 

E(P) = -tr [(tv A P) log 2 (tr A P)] (6) 

is essentailly a unique measure of entanglement ^S]- For mixed state p it seems that the basic 
measure of entanglement is the entanglement of formation |17j 

E(p) = min J2 x kE(P k ) (7) 

fc 

where the minimum is taken over all possible decompositions 

p = J2 XkPk ( 8 ) 

Again, in the case of 4 - level system, E(p) can be explicitely computed and it turns out that E(p) 
is the function of another useful quantity C(p) called concurrence, which also can be taken as a 
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measure of entanglement |11II12| . Since in the paper we use concurrence to quantify entanglement, 
now we discuss its definition. Let 

p f = (er 2 <E> cr 2 ) p(<7 2 <8> cr 2 ) (9) 
where ~p is the complex conjugation of the matrix p. Define also 

p = (p 1 / 2 pV 1/2 ) 1/2 (10) 

Then the concurrence C (p) is given by ^2 El 

C(p) = max (0,2p roax (p) - trp) (fl) 

where p m ax(/5) denotes the maximal eigenvalue of p. The value of the number C(p) varies from 
for separable states, to 1 for maximally entangled pure states. 



3 Decay in a system of closely separated atoms 

We study the spontaneous emission of two atoms separated by a distance R small compared to 
the radiation wavelength . At such distances there is a substantial probability that the photon 
emitted by one atom will be absorbed by the other. Thus the dynamics of the system is given by 
the master equation ^2] 

^ = Lp, P e £ AB (12) 

with the following generator L 

Lp =^ [2o A pa A + 2o B p o% - [p A o A + a* a B ) p - p (a A o A + o% a B )} + 

(13) 

| [2a A pa B + 2a B p a A - (a A a B + a B a A ) p - p (a A a B + af a A )] 

where 

er± =o-±®I, cr| =I®<r±, er± = ~{a 1 ±ia 2 ) (14) 

Here 70 is the single atom spontaneous emission rate, and 7 = 570 is a relaxation constant of 
photon exchange. In the model, g is the function of the distance R between atoms and g — > 1 
when R — > 0. In this section we investigate the time evolution of the initial density matrix p of the 
compound system, governed by the semi - group {T t } t >o generated by L. In particular, we will 
study the time development of entanglement of p, measured by concurrence. 

Assume that the distance between atoms is so small that the exchange rate 7 is close to 70 and 
we can use the approximation g = 1. Under this condition we study evolution of the system and 
in particular we consider asymptotic states. Direct calculations show that the semi - group {T t } 
generated by L with g — 1 is relaxing but not uniquely relaxing i.e. there are as many stationary 
states as there are initial conditions. More precisely, for a given initial state p = (pjk), the state 
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p{t) at time t has the following matrix elements 



e^pn 



Pi 2 (t) = g [^ 270t (-°i2 + Pis) + e- 70 *(pi2 - pis)] 



P22(i) 



-270* 



2 [£ 
e- 70 Vi4 



(pi2+Pl3) + e 70 *(/5l3-Pl2 



= -re 



-2 7 ot/ 



(p 22 + P33 + 2Rep 23 ) + 70 *(p 22 - P33) + 7ote 27ot pn + 7(^22 + P33 - 2Rep 23 ) 



p 23 (*) = \e- 2 ^\p 22 + p 33 + 2Rep 23 ) + ^e- 70t (p 23 - p 32 ) + l^ 2lat Pn - \{p22 + P33 - 2Rcp 23 ) 
P2i{t) = -e- 2 ^ ot (pi2 + P13) + ie-'' *(2p 12 + 2p 13 + p 24 + p 34 ) + l(p 24 - p 34 ) 

P3 3 (t) = \er 2 ^\p 22 + p 33 + 2Rep 23 ) - ^e-' at (p 22 - p 33 ) + lo t e - 2 '< at p n + l -{p 22 + p 33 - 2Rcp 23 ) 



Pu(t) 
p 44 (t) 



-270 1 



1 



(P12 + P13) + 70 *(2pi 2 + 2p 13 + p 24 + p 34 ) - -(p 24 - p 34 ) 



.I e -27ot (1 + pn _ pAA + 2 Rep 23 ) - 2 7o te- 27ot pii + 1(1 + p n + p 44 + 2Rep 23 ) 



and remaining matrix elements can be obtained by hermiticity condition pkj = Pjk- I n the limit 
f^oowc obtain asymptotic (stationary) states parametrized as follows 



where 



/o \ 

a —a /3 

—a a —0 

\0 -0 l-2aj 



a = 1(P22 + P33 - 2Rcp 23 ), = l(p 24 - p 34 ) 



(15) 



(16) 



We can also compute concurrence of the asymptotic state and the result is: 

Concurrence of asymptotic state of the semi - group {T t } generated by L with g = 1 equals to 



C(p as ) = 2\a\ = -|p 22 +p 33 - 2Rep 23 | 



(17) 



where pjk are the matrix elements of the initial state. 



4 Some examples 

In this section we consider examples of initial states and its evolution. 

I. Pure separable states 

Let 

P = = ® (18) 

where 
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are normalized. Then one can check that 



a = i (i-|(*,$)| 2 ) 

/ 3=I(|$ 2 |2^ 1 ^ 2 _|^ 2 |2 $1 $ 2 ) 



(19) 



where ( -, • ) is the inner product in C 2 . So 



C( Pas ) = -(l-|<*,<S>| 2 ) 



(20) 



From the formula (21) we sec that there are separable initial states for which asymptotic states 
are entangled. In particular, the asymptotic state has a maximal concurrence if vectors ^ and $ 
are orthogonal and its concurrence is zero (the state remains separable) if |( vp, $ ) | = 1. 
Now we discuss some special cases. 

a. When one atom is in excited state and the other is in ground state 



* = |1), $ = |0) 



the asymptotic (mixed) state is given by 



(0 








0\ 





1 

4 


1 

4 








1 
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1 

4 
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It can also be shown that in this case the relaxation process produces the state p t with concurrence 



C(pt) 



1 



-lot 



increasing to the maximal value equal to 1/2. Thus two atoms initially in separable state become 
entangled for all t and the asymptotic (steady) state attains the maximal amount of entanglement, 
b. When two atoms are in excited states 



the asymptotic state equals to 



* = $ = |1) 



|0> ® |0) 



Thus the relaxation process brings two atoms into ground states, 
c. The state |0) <g> |0) is stationary state for semi - group {T t }. 



II. Pure maximally entangled states 

Let 



P = Q(a,e 1 ,e 2 ) = 



a\/l — a' 2 i9-i 



aVl—a 2 i6 



_a 2 _ e i(6 1 +d 2 ) 
. 2 



aVl — a 2 — i9i 
2 e 



hzsl e -i(ei-e 2 ) 



aVl—a 2 c -i8?. 



1-a 2 



aVl — a 2 i8i 



i+e 2 )\ 



aVl-a 2 -i8 2 
2 e 



aVl — a' 2 „ — i0i 
2 e 
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where a £ [0,1], 61,62 € [0,27r]. Pure states Q(a, 61,62) are maximally entangled and form a 
family of all maximally entangled states of the 4 - level system 19 . It turns out that p as is defined 

by 

a=-(l-a 2 )(l~cos{6 l -6 2 )) 

\ , (21) 

and 

C(p M ) = ~(1 - a 2 )(l - cos(0! - 2 )) (22) 

From the formula (23) we see that there are initial maximally entangled states which asymptotically 
become separable (a = 1 or 6\ — 6 2 = 2kn) and such that the asymptotic concurrence is greater 
then 0. States with a = and 6\ — 6 2 = (2k + 1)tt remain maximally entangled. For example the 
state 

-L(|0>®|1)-|1)®|0» (23) 
is stable. On the other hand, the concurrence of 

-^(|0>®|1> + |1>®|0» (24) 

goes to zero faster then the concurrence of 

-L(|0)®|0> + |1>®|1» (25) 

as shown on Fig. 1. below. In Dicke's theory of spontaneous radiation processes the state (24) 
is called subradiant whereas the state (25) has half the lifetime of a single atom and therefore is 
called superradiant p^. We see that the time-dependence of concurrence reflects the relaxation 
properties of those states. 
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Fig. 1. Concurrence as the function of time for initial states: 775GO) ® |0) + |1) C§) |1)) 
(dotted line) and ^(|0) ® |1) + |1) ® |0» (solid line). 



III. Some classes of mixed states 

a. Bell - diagonal states. Let 

PB = Pl |$+)(*+| +p 2 |*-)(*-| +p 3 |* + )<* + | +P4|* _ X*-| 
where Bell states ^ and , J' ± are given by 

<&± = -L(|o) ® |o) ± |i) ® |i)), *± = ® 10) ± |o) ® |i)) 



(26) 



(27) 



It is known that all pi £ [0,1/2], pe is separable, while for pi > 1/2, pe is entangled with 
concurrence equal to 2p\ — 1 (similarly for P2, P3P4) j2(J| . Now the asymptotic state has the form 



(28) 



with concurrence C(p as ) = p^. So even when the initial state is separable, the asymptotic state 
becomes entangled. 



/o 








\ 





Pi 
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Pi 

2 








Pi 
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Pi 

2 





v° 








I -Pa) 



b. Werner states [2]. Let 



Pw = (1-p)j+p\^ + )(^ + \ 



(29) 
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If p > 1/3, pw is entangled with concurrence equal to (3p — l)/2. On the other hand 



o 







1— p p— 1 

p-1 

8 

ii 







8 u 



3+£ 

4 



(30) 



has the concurrence C(p as ) = ij 2 , so the asymptotic states are entangled for all p ^ 1. Notice 
that even completely mixed state x evolves to entangled asymptotic state. 

c. Maximally entangled mixed states . The states 



pu 



(h{5) S/2\ 

1 - 2/i(c5) 



\8/2 h(8)J 



h(S) 



1/3 6 G [0,2/3] 
(5/2 5 G [2/3,1] 



(31) 



are conjectured to be maximally entangled for a given degree of inpurity measured by tr p 2 
According to (18) the concurrence of the asymptotic state is given by 



C(Pa 



1 



(l-2h(8)) 



(32) 




0.4 



0.3 0.4 0.5 0.6 0.7 0.8 0.9 



1 trp^t 



Fig. 2. Concurrence of pm (solid line) and p as (dotted line) as the function of trpj^ 



Even in that case, there are initial states (for sufficiently small trp^j) such that the asymptotic 
state is more entangled (see Fig. 2.). 



5 Remarks on general case 



In the case of arbitrary distance between the atoms i.e. when g e [0,1), semi-group generated 
by L is uniquely relaxing, with the asymptotic state |0) ® |0). Thus, for any initial state p, the 
concurrence C{pt) approaches when t — > oo. But it docs not mean that the function t — > C(pt) is 
always monotonic. The general form of C(p t ) is rather involved, so we consider only some special 
cases. 

1. Let the initial state of the compound system be equal to |0) <£> |1). This states evolves to 



Pt 



/0 

ie-f"* (cosh 7 i + 1) 

-ie" 7ot sinh7t 

\0 



2 e 





\eT lat sinh jt 
(cosh 7 t — 1) 




(33) 



1 - e 70 * cosily 



with concurrence 



In the interval [0,i 7 ], where 



C(p t ) = e~ 7ot sinh7f 



. 1 , 7o 
t 7 = — In 

27 7o 



•7 



7 



the function (34) is increasing to its maximal value 

Wo+7 



(34) 



Cm ax 



7 



7o-7 \7o - 7 



70+7 
^7 



whereas for t > t 7 , C(pt) decreases to 0. Thus for any nonzero photon exchange rate 7, dynamics 
given by the semi - group {T t } produces some amount of entanglement between two atoms which 
are initially in the ground state and excited state. Note that the maximal value of C(pt) depends 
only on emission rates 70 and 7. 



2. For the initial states 



* ± = - 7 =(|0)®|l)±|l)®|0)) 



the relaxation to the asymptotic state |0) ® |0) is given by density matrices 



/o 





Pt 



I p -(lo±j)t 



l e -(7o± 7 )t) 





-(70±7)t 



I p -(7o±7)t 



(35) 



1 _ g-(7o±7)*^ 



with the corresponding concurrence 



C(pf) 



"(7o±7)* 



The state * is no longer stable (as in the case of 7 = 70), but during the evolution its concurrence 
goes to zero slower than C(pf) (Fig. 3.). For 7 close to 70, is almost stable. 
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Fig. 3. C(pt) (dotted line) and C(pJ) (solid line) for 7/70 = 0.99 
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